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Abstract. 

The BeppoSAX archive is currently the largest reservoir of high sensitivity simultaneous soft and hard-X ray 
data of Seyfert galaxies. From this database all the Seyfert galaxies (105 objects of which 43 are type I and 62 
are type II) with redshift lower than 0.1 have been selected and analyzed in a homogeneous way (Dadina 2007). 
Taking advantage of the broad-band coverage of the BeppoSAX MEGS and PDS instruments (~2-100 keV), the 
X-ray data so collected allow to infer the average spectral properties of nearby Seyfert galaxies included in the 
original sample and, most notably: the photon index (r~1.8), the high-energy cut-off (Ec~290 keV), and the 
relative amount of reflection (R~1.0). The data collected have been also used to test some basic assumptions of 
the unified scheme for the active galactic nuclei. The distributions of the isotropic indicators used here (photon 
index, relative amount of reflection, high-energy cut-off and narrow FeKa energy centroid) are similar in type I 
and type II objects while the absorbing column and the iron line equivalent width significantly differ between the 
two classes of active galactic nuclei with type II objects displaying larger columns (N£r~3.7x 10^^ and 6.1x10^^ 
cm~^ for type I and II objects respectively) and equivalent width (EW~220 and 690 eV for type I and II sources 
respectively). Gonfirming previous results, the narrow FeKa line is consistent, in Seyfert 2, with being produced 
in the same matter responsible for the observed obscuration. These results, thus, support the basic picture of 
the unified model. Moreover, the presence of a X-ray Baldwin effect in Seyfert 1 has been here measured using 
the 20-100 keV luminosity (EW(xL(20-100)"°'^^='=° °^). Finally, the possible presence of a correlation between the 
photon index and the amount of reflection is confirmed thus indicating thermal Comptonization as the most likely 
origin of the high energy emission for the active galactic nuclei included in the original sample. 
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1. Introduction 

The high energy emission from active galactic nuclei 
(AGN) is thought to come from the innermost regions of 
accretting systems that are centered around super-massive 
black- holes (SMBH) . For this reason, X-rays are expected 
to be tracers of the physical conditions experienced by 
matter before disappearing into SMBH. Moreover, thanks 
to their high penetrating power, energetic photons, es- 
caping from the nuclear zones, test the matter located 
between their source and the observer. Thus, they offer 
powerful diagnostics to understand the geometry and the 
physical conditions of the matter surrounding the SMBH. 

The broad-band of BeppoSAX offered for the very first 
time the opportunity to measure with a remarkable sen- 
sitivity, the spectral shape of AGN in the 0.1-200 keV 
range. This potential has been previously exploited to 
study in details a number of sample selected in different 
manners (see for example Maiolino et al. 1998; Perola et al. 



2002). These studies were fundamental in making impor- 
tant steps forward in the comprehension of the emitting 
mechanism at work in the production of X-rays (Perola et 
al. 2002) and to partially unveil the geometry of the cold 
matter surroundings the central SMBH (Maiolino et al. 
1998; Bassani et al. 1999; Risaliti et al. 1999). 

The BeppoSAX database full potential, however, was 
never exploited before. In a previous paper, the entire cat- 
alog of the Seyfert galaxies at z<0.1 contained in the 
BeppoSAX archive has been presented (Dadina 2007). 
This sample was selected starting from the catalog of 
Seyfert galaxies contained in the Veron-Cetty & Veron 
(2006) sample of AGN and contains 13 radio-loud objects 
and 7 narrow- line Seyfert 1. 

The spectral analysis was performed in the 2-100 keV 
band whenever possible and the data were fit with a set of 
template models to obtain a homogeneous dataset. Here 
the X-ray data so collected are statistically inspected in 
order to infer what are the average characteristics of the 
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nearby Seyfert galaxies contained in this sample in the 
2-100 keV band. Finally, present dataset is used to per- 
form simple tests on the imificd scheme (UM) for the AGN 
(Antonucci 1993) and on the emission mechanism acting 
in the core of the Seyfert galaxies. More detailed anal- 
ysis on this latter topic will be presented in another pa- 
per (Petrucci, Dadina & Landi, submitted) where detailed 
thermal Comptonization models (Poutanen & Svensson 
1996, Haardt & Maraschi 1993 ) will be used to fit the 
BeppoSAX data with the main scope to study the depen- 
dence of the spectral properties in the "two phase" sce- 
nario (Haardt 1991, Haardt & Maraschi 1993) assuming 
different geometries of the corona. 

2. Mean X-ray Properties 

Scope of this section is to determine and study the mean 
X-ray spectral properties of the sample and to use these 
values to test the UM model for AGN (Antonucci 1993). 
The key parameters are the ones that describe the con- 
tinuum and the absorption properties. In the framework 
of the UM models for AGN, the continuum shape is ex- 
pected to be independent from the orientation angle under 
which the source is observed. Thus, no difference should 
be measured in the parameters describing the continuum 
between type I and II objects. On the contrary, the absorb- 
ing column intervening in Seyfert 2 should be the princi- 
pal discriminator between the two classes of objects. Thus 
measuring the mean X-ray properties means to test the 
basic assumptions of the UM. 

2.1. Methods 

The origins of the X-ray photons from AGN are thought to 
be due to Comptonization of optical-UV radiation, coming 
from the accretion disk and Comptonized by the e~ in the 
hot corona that sandwiches the disk (Haardt & Maraschi 
1991; Haardt 1993, Poutanen & Svensson 1996, Czerny 
et al. 2003). The mechanism is assumed to be. at least 
at the zcro-th order, very similar in each Seyfert galaxy. 
Under this hypothesis, the differences between the X-ray 
spectra of different objects are supposed to be mainly due 
to two kind of factors: i) the time-dependent state of the 
emitting source; ii) the intervening matter that imprints 
on the emerging spectrum the featTircs typical of its phys- 
ical state. In such a scenario, the observations of many 
sources can be regarded as the long-term monitoring of a 
single source. On the other hand, it is also true that the 
contrary has some comparison in the literature: e.g. time 
sparse observations of a single source in different states 
can resemble observations of objects with completely dif- 
ferent spectral properties. This is the case, for example, of 
the narrow line Seyfert 1 NGC 4051 that displayed vari- 
ations in flux/luminosity by a factor of ~100 (Guainazzi 
et al 1998) associated to strong variations of the spec- 
tral shape (r~0.5-2.4; Guainazzi et al. 1998; Turner et al. 
1999; Ponti et al. 2006, but see also Crenshaw & Kraemer 
(2007) for a different interpretation of this spectral behav- 



ior in terms of variable ionized absorption). To calculate 
the mean X-ray properties of the sources included in this 
work, thus, I treated the multiple observations of single 
sources separately: i.e. I assumed different observations of 
the same source as if they were observations of different 
sources. 

This method, in principle, is expected to work prop- 
erly for all those quantities which are supposed to vary 
in accordance with the state of activity of the central 
nucleus. For example, the photon index F is known to 
vary accordingly with the AGN flux state (Lee et al. 2000; 
Shih, Iwasawa & Fabian 2002, Ponti et al. 2006) and the 
high-energy cut-off (Ec) is linked with the temperature 
of the corona and thus expected to be variable (Haardt, 
Maraschi & Ghisellini 1997). On the contrary this method 
is not expected to work properly when constant compo- 
nents are considered. This could be the case, for instance, 
of the cold absorption assumed to be due to the putative 
dusty torus (Antonucci 1993). For this component, thus, 
the average value recorded for each source should be used. 
Nonetheless, also these components were observed to vary 
in a number of objects (see for example the case of NGC 
4151, De Rosa et al. 2006 or Risaliti et al. 2002). Moreover, 
the constancy of the properties of the cold absorption is 
predicted under the hypothesis of a continuous distribu- 
tion of the matter that forms the torus. On the contrary, if 
the torus is formed by blobs/clouds (Elitzur & Shlosman 
2006), variability in the measurements of absorbing col- 
umn is naturally expected. For both these observational 
and theoretical reasons also the N^f measured in each sin- 
gle observation were treated separately. 

In between these two cases, arc the properties of the 
emission FeKa line. BeppoSAX had a too low sensitiv- 
ity to detect the relativistically broadened component of 
this feature in a large number of sources. Thus it turned- 
out that only the narrow line has been detected in the 
vast majority of the objects included in the original cata- 
log. This component is supposed to originate far from the 
SMBH, at least at --1000 Schwarzschild radii (Mattson & 
Weaver 2004), i.e. very likely at the inner edge of the torus 
(Nandra 2006). At these distances from the SMBH, the 
relativistic effects are negligible. Thus, in principle, the 
parameters that describe the line should be regarded as 
stable. Nonetheless also the properties of the narrow FeKa 
line were observed to vary with its line energy centroid 
changes between ^^6.4 keV (neutral Iron) and ^6.9 keV 
(H-like Iron) showing different ionization states. Moreover, 
the equivalent width (EW) of this component is not only 
function of the intensity of the line itself, but it is also di- 
rectly linked with the underlying variable continuum emit- 
ted in the regions close to the SMBH. This makes EW a 
variable quantity. The parameters of the emission FeKa 
line were, thus, treated as variable ones in order to test 
its origin, and not averaged when objects were observed 
more than once. 

Finally, it is worth noting here that in a number of 
cases, it was necessary to deal with censored data (see ta- 
ble 1). To manage these data properly, the ASURV soft- 
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ware (Feigelson & Nelson 1985; Isobe et al. 1986) has 
been used. In particular, to establish if the distributions 
of parameters of type I and type II objects were drawn 
from different parent populations, the Peto & Prentice 
generalized Wilcoxon test (Feigelson & Nelson 1985) has 
been used while to calculate the mean values consider- 
ing also the censored data the Kaplan-Meyer estimator 
has been used. To establish the presence of correlations 
between different quantities, both the Spearman's p and 
the generalized Kendall's r methods were applied. The 
linear regressions were calculated using the Bukley-James 
and Schmitt's methods. In the following, two quantities 
are considered as drawn from different parent population 
when the probability of false rejection of the null hypoth- 
esis (same parent population) is Pn«H<l%- Similarly, one 
accepts that there is a correlation between two given quan- 
tities when the probability of absence of correlation re- 
mains smaller than 1%. 

Table 1. General characteristics of the data analyzed in 
this work. The number of detections and censored data are 
reported for the interesting parameters for the whole sam- 
ple of objects (columns 2 and 3), for the Seyfert 1 galaxies 
(columns 4 and 5), and for the Seyfert 2 objects (column 
6 and 7). The 90% confidence interval limits were used 
for censored data and the detected values were defined if 
determined with a 99% confidence level (Dadina 2007) 



Parameter 


Tot. 


Sample 


Seyfert 1 


Seyfert 2 




Det. 


Cens. 


Det. 


Cens. 


Det. 


Cens. 


Nff 


83 


80 


31 


53 


52 


27 


EW 


129 


7 


66 


4 


63 


3 


R 


68 


18 


46 


9 


22 


9 


Ec 


33 


51 


27 


26 


6 


25 



2.2. The X-ray continuum and the cold absorption 

The X-ray continua of the sources have been modeled us- 
ing a cut-off power-law, that describes the primary emis- 
sion from the hot corona plus a reflection component 
{PEXRAV model in Xspec, Magdziarz & Zdziarski 1995) 
to account for the contribution expected to be due to the 
disk. It is worth noting, however, that an additional re- 
flection component could rise from the torus (Ghisellini, 
Haardt & Matt 1994) and that the disentangling between 
the two reflection component is impossible with the qual- 
ity of the available data. Wherever its origin, the reflection 
has been assumed to be due to cold matter. The inter- 
esting parameters are the photon index (F), the relative 
amount of reflection (R), and the high-energy cut-off (Ec). 



In table 2 the results for the whole set of observations 
and for the two classes of Seyferts are reported as well as 
the probability that Seyfert 1 and 2 are drawn from the 
same parent populations. The histograms of the distribu- 
tions of the interesting parameters for the entire sample 
of objects (first column) and for type I (second column) 
and type II (third column) are reported in figure 1. 

Table 2. Mean spectral properties. Col. I: Spectral pa- 
rameter; Col. II: Seyfert 1 mean value; Col. Ill: Seyfert 2 
mean value; Col. IV: Probability that Seyfert 1 and Seyfert 
2 are drawn from the same parent populations. 



Parameter 


Tot. 


Seyfert 1 


Seyfert 2 




F 


1.84±0.03 


1.89±0.03 


1.80±0.05 


90% 


R 


1.01±0.09 


1.23±0.11 


0.87±0.14 


5% 


Ect 


287±24 


230±22 


376±42 


5% 




31.7±9.1 


3.66±2.34 


61.3±18.0 


<0.1% 



^ in units of eV; * in units of 10^^ cm 



As previously said, the UM for AGN (Antonucci 1993) 
predicts that F, R, and Ec are observables independent 
from the inclination angle, thus the two classes of Seyfert 
galaxies should display very similar characteristics. This is 
confirmed by the analysis of the present sample. In partic- 
ular, there are no hints that the distributions of photon in- 
dex F for the two types of Seyfert are drawn from different 
parent populations (P„u(i^90%). Moreover, the photon- 
index peaks, for both classes, between 1.8-1.9 well in 
agreement with the two-phase models for the production 
of the X-ray in Seyfert galaxies that predicts F^l.5-2.5 
(Haardt & Maraschi 1991; Haardt 1993: Haardt, Maraschi 
& Ghisellini 1997). Few objects have extremely flat spectra 
with F<1. Type II objects that show so hard X-ray spec- 
tra are supposed to be Compton-thick sources for which, 
in the 2-10 keV band, only the reflected/flat spectrum is 
observable. This is the case for NGC 2273 which displays 
the harder X-ray spectrum. This source was first classi- 
fied as a Compton-thick object by Maiolino et al. (1998). 
The Seyfert Is with fiattest spectra are NGC 4151 that is 
known to have a hard spectrum with complex and variable 
absorption (De Rosa et al. 2006) and Mrk 231 (classified 
as a type I AGN by Farrah et al. 2003). The latter source 
shows a very hard X-ray spectrum with F~0.7. This source 
is also classified as a BAL QSO (Smith et al. 1995). A re- 
cent spectral analysis in the X-ray band of this source was 
presented in Braito et al. (2004). By combining XMM- 
Newton and BeppoSAX data, these authors speculated 
that the spectrum of the source below ~10 keV is refiec- 
tion dominated, thus presenting a case that can be hardly 
reconciled with the UM of AGN. Moreover, Mrk 231 is 
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Log(Nh) Log(Nh) Log(Nln) 



Fig. 1. Photon index F {first row), R (second row), Ec [third row), and N/j (fourth row) in units of cm ^ distributions for the 
whole dataset (left column), for type I objects (center column), and for type II objects (right column). 



an ultra-luminous infra-red galaxy. These sources display 

strong starburst activity that can dominate the total X- 
ray luminosity of the galaxy (see Franceschini et al. 2003 
and Ptak et al. 2003 for details on this topic) although 
in Mrk 231 at least 60% of the observed 0.5-10 keV flux 
seems to be due to the AGN component (Braito et al. 
2004). 

Less conclusive results arc obtained for R and Ec. In 
both cases, the probability of false rejection of the null hy- 
pothesis (the two distributions are drawn from same par- 
ent populations, in accordance with the UM predictions), 
is PrmH~5%. It is worth noting, however, that both type I 
and II objects lay in the same range of R and Ec. In partic- 
ular, the distribution of R for type I objects is dominated 
by a peak of detections between R~l-1.2 due to the contri- 
bution of a few sources observed several times (e.g. NGC 
5548, IC 4329a). This has probably introduced a system- 
atic effect not smeared-out by the relatively small number 
of useful observations (i.e, the ones for which R and Ec 
have been estimated). On the other hand, the result on 
Ec is most probably polluted by the large number of lower 
limits for Seyfert 2. In fact if the same test is performed 
only on the detections it is obtained that PnMZi~35%, mak- 
ing impossible to assess that the two samples are drawn 
from two different populations. 



The result is unambiguously in accordance with UM 
predictions for the absorbing column. In this case the sta- 
tistical tests confirm that Seyfert 1 and Seyfert 2 dis- 
play very different absorption characteristics with the type 
II objects being more heavily absorbed than Seyfert 1. 
Again, it is noticeable that a few type I objects show high 
column densities (up to ~10^^ cm~^) and some Seyfert 2 
have low columns (down to lO^'^ cm~^). These are not new 
results: the high column of NGC 4151 (~5xl022, De Rosa 
et al. 2006) is well known. The highest column measured 
in a Seyfert 1 is detected during the Jime 9, 1998 obser- 
vation of NGC 4051. During this observation the source 
appeared "switched-off " and only a pure reflection compo- 
nent was measured (Guainazzi et al. 1998). The spectral 
fit in Dadina (2007) degenerated between two-solutions, 
one in which the source was purely reflection-dominated 
(R>7) and a second one in which a direct component 
was visible but highly absorbed. The latter scenario was 
slightly preferred by a pure statistical point of view when 
the 2-200 keV band is condidered, and, for homogeneity, 
entered in the catalog (Dadina 2007). Nonetheless, when 
the entire BeppoSAX band (0.1-200 keV) is considered, 
the reflection scheme is preferred (Guainazzi et al. 1998). 

Finally, a number of objects show upper-limits of the 
order of ^10^^ cm~^ to the absorbing column. This is a 



Mauro Dadina: Average X-ray properties of Seyferts observed by BeppoSAX 



5 




Fig. 2. Distribution of the Feka emission line energy centroid {first row) and of its EW (second row). Distributions for the whole 
sample of observations (left panel), for the type I objects (center panel) and for type II objects (right panel). 



selection effect induced by the energy band considered us- 
ing only MECS and PDS data (~2-100 keV, Dadina 2007). 
The low-energy cut-off due to such a column (10^^ cm^^) 
peaks in fact at E~2 keV and only for those objects with 
a good statistics it is possible to infer upper limits on the 
N/f below 10^^ cm^^. This is most probably responsible of 
the high value obtain for the average Njj in type I objects. 

3. Probing the origin of the FeKa emission line 

The FeKa line is produced by reprocessing the primary X- 
ray emission in matter surrounding the source of hard pho- 
tons. In the framework of the UM, the origin of this com- 
ponent can be placed in a number of region s such as the 
accretion disk, the dusty torus, and the broad-line regions 
(even if this last hypothesis is disfavored by the recent 
results obtained with the XMM-Newton and Chandra 
observatories and presented in Nandra 2006). If the line 
originates in the disk close to the SMBH, relativistic ef- 
fects that broaden the resulting line are expected. For the 
vast majority of the sources included in the original sam- 
ple, only narrow component of such feature were detected 
and only in a few cases broad emission lines (e.g. IC 4329a) 
or relativistically blurred features were detected (for ex- 
ample in MCG-6-30-15). Thus, the results presented here 
are essentially based on the measured properties of the 
narrow features. 

As shown in figure 2 (first row), the line energy cen- 
troid is peaked at ~6.4-6.5 keV (see also table 3) in both 
type I and II objects. The centroid is slightly above 6.4 
keV but, considering the energy resolution of the instru- 
ment at these energies (~200 eV FWHM, Boella et al. 
1997), the results obtained here are in agreement with the 
line being mainly produced in cold or nearly cold matter 
(ionization state below FeXVII), i.e. in both type I and II 
objects, by matter in the same physical state. However, a 
well known difference between the two classes of objects is 
the EW of the narrow FeKa line in type II objects, which 
shows stronger features than type I objects (see second 
row of figure 2 and table 3, Bassani et al 1999, Risaliti 



et al. 1999, Cappi et al. 2006 and Panessa et al. 2006). 
As shown in the central panel of figure 2, the Seyfert 1 
peak at EWFeifa"^100-200 eV while the type II have a 
broader distribution with an hard tail that reaches values 
well above 1 keV. Also few Seyfert 1 have large values of 
EW of the FeKa line (above; 300-400 eV). The large EW 
tail of the Seyfert 1 distribution is composed mainly of ob- 
jects in which the broad components of the FeKa line are 
detected such as in MCG-6-30-15, Mrk 841 and IC 4329a. 
The Seyfert 1 with largest EW is NGC 4051 during the 
June 9, 1998 observation when its spectrum was due to 
pure reflection (Guainazzi et al. 1998). 

The larger FeKa EW in Seyfert 2 galaxies is in agree- 
ment with the UM (Antonucci 1993). If the origin of this 
component is indeed located in the dusty torus, than the 
line EW has to be correlated with the absorber column 
density. This is indec^d what is observed also in this sample 
(see figure 3, left panel). Moreover, the Spearman p and 
Kendall's r tests indicate that a correlation between the 
FeKa EW and the N/f is highly probable for type II ob- 
jects (Pntiii<0.1%), i.e. for that sources for which we can 
have direct evidences of the torus absorbing column. The 
robustness of the Njj estimates have been tested by corre- 
lating it with the model independent indicator offered by 
the ratio of the observed fluxes at 2-10 and 20-100 keV re- 
spectively (center panel of figure 3). The two quantities are 
strongly correlated (generalized Spearman p and Kendal 
T tests give Pntiii^O.1%) with only a few exceptions: the 
pure Compton-thick sources which are located in the di- 
agram below the majority of the sources. This effect is 
expected since the absorbing column can affect the X-ray 
radiation below ~10 keV while at harder energies the ra- 
diation pierces the matter for columns <3-5xl0^^ cm~^. 
Otherwise, the Compton absorption dominates and also 
photons with energy above 10 keV are stopped since the 
Klein-Nishina regime is reached. 




Fig. 3. Left panel: Log(EWFe/fct) vs. Log(N£f). As expected, the sources are divided in two families: the ones that follow the 
expectations if the FeKa line is produced in the absorption matter and the candidate Compton-thick ones that display low 
absorption and large EW. The solid line indicates the prediction by Makishima (1986). Center panel: Log(F2-iofeev/F2o-ioofceK) 
vs. Log(N/f ). The two quantities are correlated as expected since the 2-10 keV band is strongly affected by the absorption while 
the 20-100 band is almost free from absorption. The solid line is the best fit obtained with linear regression methods Right Panel: 
Log(EWFeKa) VS. Log(F2^ lofcev /F20- loofee V ) • The two quantities are strongly correlated ( Pnu((<0.01). The ratio between X-ray 
fluxes is a model independent indicator of the absorption affecting the low energy band. Thus, this correlation strongly indicates 
that the narrow component of FeKa line in emission is indeed produced by the same matter responsible for the absorption. The 
solid line is the linear regression obtained using Bukley- James method (Isobe et al. 1986). 



Table 3. Mean properties of the FeKa emission line. Col. 
I: spectral parameter; Col. II mean value for the whole 
sample; Col. Ill: mean value for Seyfert 1; Col. IV: mean 
value for Seyfert 2. Col. V: Probability that the parame- 
ters of type I and type II objects are drawn from the same 
parent population. 



Parameter 


Tot. 


Seyfert 1+ 


Seyfert 2+ 






6.49±0.02 


6.46±0.03 


6.51±0.03 


32% 


EWf^kJ 


448±67 


222±33 


693±195 


<1% 



in units of keV; * in units of eV 



As stated above, when the EW of the FcKq omission 
line is tested against the measured Nh (left panel of Figure 
3) a result in good agreement with what predicted by the- 
oretical models is obtained (Makishima 1986; Leahy & 
Creighton 1993). The majority of the sources, in fact, be- 
have as expected if the line is produced by the absorb- 
ing matter that depress the direct continuum (Makishima 
1986). All the known Compton-Thick sources are located 
at low Njf and high EW, in accordance with previous re- 
sults (Bassani et al. 1999, Risaliti et al. 1999 ). As an 
additional test, the EW of the FeKa line has been plotted 
versus the F2_iofcey/F2o-ioofceV ratio. As expected (see 
right panel of figure 3), a good correlation ( Pn«n<0.1% 
according to generalized Spearman p and Kendall r tests) 
is obtained. 



These results thus confirm that the properties of the 
FeKa line agree with the expectations of the UM for 
AGN. Nonetheless, this is not the only information we 
have about the Iron line. In recent papers (Iwasawa & 
Taniguchi 1993; Page et al. 2004, Grandi et al. 2006; 
Bianchi et al. 2007) it has been claimed that a X-ray 
"Baldwin effect" (or Iwasawa- Taniguchi effect) is present 
in AGN when the FeKa intensity is probed against the 
2-10 keV luminosity. Here this effect is tested considering 
for the first time both the 2-10 keV and the 20-100 keV 
luminosities. 

A strong correlation ( Pnu/i^O.l'X; using Spearman p 
and generalized Kendal r tost) is found when the EW of 
the FeKa line is plotted both against the observed 2-10 
keV (panel (a) of figure 4) and 20-100 keV (panel (c) of fig- 
ure 4) keV luminosities. The nature of these correlations, 
however, is not straightforward, especially for the 2-10 keV 
luminosity. In this energy band the effect of the absorber 
is very important. As previously demonstrated, the EW of 
the Iron line correlates with N^, but stronger absorbing 
columns imply lower fluxes. Moreover, when the relation 
between the FeKa line EW and the 2-10 keV flux is in- 
vestigated (panel (b) of figure 4), it is found that the two 
quantities are correlated {Pnuii^^%)- This is explainable 
only in terms of instrumental effects, since the sensitiv- 
ity of the instruments to narrow features decreases with 
the source's flux. Thus, to be detected in faint objects, 
the FeKa line must be strong enough. Since the original 
sample is limited to the local Universe, this correlation in 
flux acts, at least partially, also in the EW vs L2_iofeeV 
relation. 
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38 40 42 44 46 41 42 43 44 45 46 -12 -11 -10 -9 -8 

Log(L20-100) Log(L20-100) Log(F20-100) 

Fig. 4. Panel (a): Log(EWFeKc) vs. h2-io,observed- Panel (b): Log(EWFeKo) vs. F2-io,o6sert.ed. Panel (c): Log(EWFejf „ ) vs. 
L20-100, observed. Panel (d): The same of panel (c) but only for Seyfert 1 objects. Panel (e): Log(EWFeK„ ) vs. F20-100, observed. 
Solid lines in panel (a), (b), and (c) are the linear regressions obtained for the whole sample of observations. The dashed line 
in panels (c) and (e) is the linear regression obtained considering only type I objects. As visible in panel (b) the less scattered 
relation is obtained considering the 2-10 keV observed flux. The relation is linear, as expected if the correlation is due to selection 
effects, i.e.considering that in faint objects only large EW were detectable by the MECS instruments on-board BeppoSAX. 



This Nfl- effect should be negligible when the 20-100 
keV band is considered. In fact, in this case, one expects 
to find a correlation between the observed 20-100 kcV lu- 
minosity and the EW of the FeKa line only in the most 
extreme cases, i.e. for the "pure Compton thick" objects. 
Apart from NGC 1068, these sources are too weak to be 
detected by the PDS, thus unable to drive the relation ob- 
served in plot (c) of figure 4. Moreover, panel (e) of figure 4 
indicates that the EW of the iron line is not related to the 
20-100 keV flux thus excluding that the "Baldwin effect" 
measured using the 20-100 keV band is due to instrumen- 
tal selection effect as it happens for the F2-iofeeV- On the 
other hand, in the 20-100 keV band the reflection-hump at 
~30-40 keV contributes to the observed flux. If the origin 
of the FeKa line is due to the same matter responsible of 
the reflection, one should expect that the EW of the Iron 
line should increases as the reflection component augments 
the 20-100 keV flux (i.e. the Iron line EW should correlate 
with the 20-100 keV flux). Thus, the net contribution of 
the reflection component should act in the opposite direc- 
tion to that observed (i.e. larger iron line EW at smaller 



20-100 keV flux). If the two classes of Seyfert galaxies are 
analyzed separately it is obtained that a strong correlation 
is found for Seyfert 1 (Pmtii<0.1% using both Spearman p 
and generalized Kendal t tests, panel (d) of figure 4) while 
no correlation is evident for type II objects (PnM«~80%). 
This result is not unexpected since the EW of the obscured 
sources are boosted by the suppression of part of the un- 
derlying continuum. To conclude, the presence of a X-ray 
Baldwin effect for Seyfert 1 is unambiguously confirmed 
by present data if the L2o-ioofeeV is considered and it has 
the following relation: 

Log(EW)=-(0.22±0.05)xLog(L2o-ioo)+11.91±2.52 (1) 

The slope of the relation found in this work is in agree- 
ment with what previously obtained by Page et al. (2004) 
(EWa(L2_io)~'''^^''^° °*) using a sample containing both 
radio-quiet and radio-loud objects. Present result is also 
consistent with what found by Zhou & Wang (2005) (who 
used a sample containing both radio-quiet and radio-loud 
objects founding EWoc(L2-io)-°-i^=^°-°5) and Bianchi et 
al. (2007) (who used only radio-quiet objects obtaining 



8 Mauro Dadina: Average X-ray properties of Seyferts observed by BeppoSAX 




Log{Nh} Log(Ec) Gamma 

Fig. 5. Left panel: Photon index V plotted versus the measured Nif (in units of cm^). No significance trend relating these two 

quantities is found; center panel: R vs. Ec (in units kcV), no correlation is found between these quantities; Tight panel R vs. V. 
The two quantities are correlated with a high significance level (Pn«ii<l%). The line is the linear regression best fit. 



EWoc(L2-io)°-^^=^° °^). On the contrary, Jiang, Wang & 
Wang (2006) found that, excluding the radio-loud AGN 
from a sample similar to the one used by Page et al. 
(2004), found that EWoc(L2-io)°-^°=^° °^ It is worth re- 
calling here, that the the present sample is composed by 
both radio- loud and radio quiet sources (Dadina 2007). In 
particular, it contains 7 radio-loud Seyfert 1, and only 
for 5 of them the 20-100 and iron line data are avail- 
able. Nonetheless, the presence of these sources in the 
sample does not affect the FeKa EW vs. L20-100 rela- 
tion (EWoc(L2-io)°'^^'''° °^ excluding radio-loud type I ob- 
jects). 

The origin of the X-ray "Baldwin effect" is unclear. In 
the light bending scenario (Miniutti & Fabian 2004) the 
height of the source above the accretion disk determines 
the degree of beaming along the equatorial plane of the 
high energy emission. Because of relativistic effects, the 
closer the source is to the disk, the greater will be the frac- 
tion of X-rays beamed in the equatorial plane (i.e. towards 
the disk) and correspondingly lower will be the observed 
flux. Thus, the EW of the relativistically blurred FeKa 
line produced in the inner regions of the disk and the EW 
of the narrow Iron line produced in the outer parts of the 
disk would appear enhanced in low-state sources. 

On the contrary. Page et al. (2004) speculated that this 
effect could indicate that luminous sources are surrounded 
by dusty tori with lower covering factor thus pointing to- 
wards a torus origin of most of the narrow FeKa line. 
The present work supports this view. The FeKa line EW 
of the Iron line correlates with the observed Njy as pre- 
dicted by theory (Makishima 1988; Lehay & Creighton 
1993; Ghisellini, Haardt & Matt 1994). Moreover, also the 
case of the extremely low state of NGC 4051 (Guainazzi 
et al. 1998) included in this dataset seems to point in this 
direction. The huge EW of the narrow FeKa line recorded 
in this observation, in fact, is typical of Compton-thick 
type II objects, but the line does not show evidence of rel- 



ativistic broadening due to the contribution of the inner 
orbits of the accretion disk. 

4. The r-R relation. 

In the fitting procedure some parameters may degenerate 
given the interdependence among them. This is the case, 
for example of the photon index with the column density 
for low statistics observations. The same could happen for 
the determination of R and Ec, since R introduces in the 
spectrum a bump peaked between 20-40 keV and declining 
at higher energies where the Ec may appear. To check if 
the results presented here are affected by such effects, the 
correlations between these parameters have been studied 
and the results are presented in Figure 5. 

Left panel of Figure 5 shows how, on average, the esti- 
mate of F is not affected by the simiiltaneous determina- 
tion of the absorbing column. In fact, no trend is observed 
between F and Njj. Obviously, this does not imply that 
this is true for each single source included in the origi- 
nal sample. On the other hand, this is an expected result 
since the broad band of BeppoSAX should reduce this 
spurious effect. Similar results are obtained also when the 
r vs. Ec, and R vs. Ec (center panel of Figure 5) rela- 
tions are investigated. In these cases the Spearman p and 
Kendall's r tests do not sustain the existence of a relation 
between these quantities (PnMi('~15%). All these indica- 
tions suggest that, if any, the possible degeneracy in the 
fitting procedure did not introduce strong spurious rela- 
tions between the spectral parameters. 

However a strong correlation (PnuH<0.01%) is 
recorded between F and R that are linked by the following 
relation: 

R=(4.54±1.15)xr-(7.41±4.51) (2) 

It is hard to define if this correlation is the result of 
a systematic effect or not since it is possible that these 
two quantities degenerate in the fitting procedure. Flat 
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power-law with small reflection could be described also by 
steep power-law plus strong reflection. The total absence 

of similar correlations between F and Ec and R vs. Ec 
seems to suggest that this correlation is indeed real. 

A similar relation was previously found using Ginga 
and RXTE data (Zdziarski, Lubinski, & Smith 1999; 
Gilfanov, Churazov & Revnivtsev 1999). Zdziarski, 
Lubinski, & Smith (1999) interpreted it as evidence of 
thermal Comptonization as origin of X-rays providing that 
the optical-UV seed photons were mainly produced by the 
same material responsible for the reflected component. In 
this case, in fact, the cooling rate of the hot corona is di- 
rectly linked to the power-law slope. But the cooling rate 
is also related to the angle subtended by the reflector. This 
result is also in agreement with predictions of models that 
consider mildly relativistic outflows driven by magnetic 
flares (Beloborodov 1999). More in general, Merloni et al. 
(2006) demonstrated that any geometry in which the hot. 
X-ray emitting plasma, is photon starved (i.e. geometries 
of the accretting systems in which the accretion disk is 
only partially covered by the Comptonizing plasma such 
as patchy coronae, inner ADAF plus outer disks etc.) will 
produce hard X-ray spectra, little soft thermal emission 
and weak reflection component. On the other hand, ge- 
ometry corresponding to a very large covering fraction of 
the cold phase, have strong soft emission, softer spectra 
and strong reflection fraction (CoUin-Souffrin et al. 1996). 
Thus, moving along the F vs. R relation it implies moving 
from lower to higher accretting systems. 



5. Summary and conclusions 

The average properties of Seyfert galaxies in the local 
Universe (z<0.1) as seen by BeppoSAX has been in- 
vestigated, analyzing the sample of objects presented in 
Dadina (2007). Multiple observations of single objects 
were treated independently, i.e. the multiple measure- 
ments of parameters were not averaged for statistical pur- 
poses. This method has been chosen since, in the frame- 
work of the simplest version of UM for AGN (Antonucci, 
1993) the AGN are thought to be very similar to each- 
other and only orientation/absorption effects and the 
activity-level of the targets could introduce observational 
differences between different objects. In this scenario, the 
monitoring of a single source could be reproduced by the 
observations of many sources in different states and "vice 
versa" . 

BeppoSAX offered the advantage of a useful X-ray 
broad energy band. Data studied here fall in the 2-100 
keV band for a majority of objects. This advantage has 
been used to investigate the properties of the high-energy 
continuum of Seyfert galaxies. As stated in the previous 
paper (Dadina, 2007), the basic template was a power- 
law with a high energy cut-off plus a reflection component 
(namely PEXRAV model in XSPEC). The results of this 
analysis can be summarized as follow: 



— the average slope of the power-law is 1.84±0.03 for the 
entire sample of objects. Considering the two families 
of AGN separately, it turns-out that F=1.89±0.03 for 
type I objects (including Seyfert 1, 1.2 and 1.5) while 
r=1.80±0.05 for type II objects (including Seyfert 1.8, 
1.9 and 2); 

— the average value of the relative reflected-to-direct nor- 
malization parameter R is 1.01±0.09 with a slight dif- 
ference between the two classes of Seyfert galaxies 
(R=1.23±0.11 for type 1 objects and R=0.87±0.14 for 
the type II ones); 

— the high energy cut-off was measured to be Ec=287±24 
(Ec=230±22 keV for Seyfert 1 and Ec=376±42 keV for 
Seyferts 2); 

— as expected and as known from previous works, the 
absorbing column is very different in the two classes 
of objects. On average N//~3. 66x10^^ cm~^ type I 
and Nh-6. 13x1023 cm^^ for type II AGN. The high 
mean value obtained for Seyfert 1 is caused by a se- 
lection effect induced by the energy coverage of the 
MECS-hPDS instruments (2-100 keV). 

— evidence of a X-ray Baldwin effect is found in Seyfert 
1 galaxies when the EW of the FeKa line is plotted 
against the 20-100 keV luminosity. 

— a significant correlation has been found between R and 
F. 

These results are well in agreement with the basic as- 
sumptions of the UM for AGN (Antonucci 1993). In fact, 
no differences are measured in the obscrvables that are 
supposed to be isotropic while the absorbing column seems 
to be the only discriminator between the different types of 
Seyfert galaxies. This reflects also in the properties of the 
FeKa line. No difference is measured in the line centroid 
(see table 3) between the two classes of Seyfert galaxies. 
Type II objects, however, display more intense features 
(EW=222±33 eV for Seyfert 1 and EW=693±195 eV for 
Seyfert 2). The physical origin of the X-ray "Baldwin ef- 
fect" here measured for Seyfert 1 using the 20-100 lumi- 
nosity is unclear. Both light bending (Miniutti & Fabian 
2004) and torus models (Page et al. 2004) are consistent 
with present data even if the strong relation of the FeKa 
line EW in type II objects with the absorption column in- 
dicates that the most of the narrow line component should 
be due to the torus. Finally, the measTired F-R relationship 
is consistent with thermal Comptonization models. 
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